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We study the 5=1 Heisenberg antiferromagnet on a spatially anisotropic triangular lattice 
by the numerical diagonalization method. We examine the stability of the long-range order of 
a three-sublattice structure observed in the isotropic system between the isotropic case and the 
case of isolated one-dimensional chains. It is found that the long-range-ordered ground state 
with this structure exists in the range of 0.7 ^ J2/J1 < 1, where Ji is the interaction amplitude 
along the chains and J2 is the amplitude of other interactions. 

KEYWORDS: Antiferromagnetic Heisenberg spin model, triangular lattice, frustration, spatial anisotropy, 
numerical diagonalization method, Lanczos method 



Frustrated magnets have attracted much attention 
from the viewpoint of realizing exotic quantum states 
and phase transitions that occur between such states. 
One such magnet is the triangular-lattice Heisenberg an- 
tiferromagnet. Much effort has been devoted to studies of 
the S — 1/2 model, particularly since Anderson^ pointed 
out that this model is a possible candidate for the re- 
alization of the spin liquid ground state due to mag- 
netic frustrations. Now, it is widely believed that the 
ground state of the S ^ 1/2 model has a long-range 
order (LRO) with a small three-sublattice magnetiza- 
tion.-^"^ This spin structure is also called the 120 deg 
structure from the directions of neighboring spins. How- 
ever, estimating the magnetic order of a three-sublattice 
structure (3SS) quantitatively by direct methods that 
are unbiased beyond any approximation or variational 
method is still difficult even now. Numerical diagonal- 
ization data of small finite-size (FS) clusters of up to 
36 sites of the S* = 1/2 model^"'^ were examined.^ How- 
ever, Leung and Runge^ and Bernu et al.^ reported the 
subtlety in quantitative extrapolation from their diago- 
nalization data of spin correlation functions for observing 
the order directly. 

Since the LRO occurs as a consequence of the deli- 
cate balance of the frustrated situation, the instability 
of this order is an attractive issue in the case when the 
model includes factors additional to the triangular-lattice 
Heisenberg antiferromagnet. One such factor is the spa- 
tial anisotropy of interactions in the system. Let us con- 
sider that the interaction along a specific direction (Ji) 
is different from the interaction along the other two di- 
rections (J2) among three equivalent directions in the 
isotropic case (see Fig. 1); here, we focus our attention 
on the range of < J2/ Ji < 1. In the case of J2 = 0, the 
system is reduced to isolated spin chains with an antifer- 
romagnetic interaction, in which the LRO does not oc- 
cur owing to the one-dimensionality. Therefore, the order 
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Fig. 1. Spatial anisotropic triangular lattice. 



must disappear at a point in this range of < J2/ Ji < 1. 

In the classical case of 5 = 00, the spiral ground state 
is realized in the range of < J2/J1 < 1; the behavior is 
characterized by the canting angle between neighboring 
spins. In the quantum case, however, not only the cant- 
ing angle but also sublattice magnetization in the spiral 
state can be changed due to quantum fluctuation; the sit- 
uation is different from the classical case. The = 1/2 
case has also been tackled. ^'^"^^ On one hand, there are 
various reports concerning the range of the spiral state 
depending on the method. On the other hand, it was 
shown by the functional renormalization group approach 
that the LRO in the isotropic case disappears quickly 
when the anisotropy is switched on.^^ Our understand- 
ing of this S* = 1/2 model has not reached a consensus 
so far. 

In this letter, we examine the S = \ triangular-lattice 
Heisenberg antiferromagnet with a spatial anisotropy. 
One major difference of the 5* = 1 case from the 5 = 1/2 
case is that the ground state at J2 = in the S — 1 case 
shows the Haldane gap, whereas that in the S = 1/2 
case is gapless in the spin excitation. It is a nontrivial 
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issue how the existence of the gap affects the magnet- 
ically ordered phase, which is supposed to be realized 
in the isotropic case. Only two approximations^*'' have 
been attempted to elucidate the S — 1 case. By a linked- 
cluster series expansion method,^" it was shown that the 
ratio J2/ Ji at the boundary between the magnetically or- 
dered phase of the spiral state and the disordered phase 
is 0.33. In ref. 21, it was shown by the coupled cluster 
method that the ordered phase spreads over a range of as 
low as J2/ Ji — 0.25. In order to obtain an understanding 
of what beyond any approximation, direct numerical sim- 
ulations have become increasingly important. However, 
the quantum Monte Carlo method cannot be applied to 
the analysis of the triangular-lattice antiferromagnet ow- 
ing to the so-called negative-sign problem from the frus- 
tration in the system. Although the density matrix renor- 
malization method is powerful for the analysis of one- 
dimensional systems irrespective of whether or not frus- 
tration exists, this method is much less effective for the 
analysis of two-dimensional systems such as a triangular- 
lattice antiferromagnet. Under these circumstances, the 
numerically exact diagonalization method based on the 
Lanczos algorithm is applicable to the triangular-lattice 
antiferromagnet irrespective of the spatial dimensional- 
ity and the existence of magnetic frustration. The pur- 
pose of the present study is to extract information on 
the magnetic LRO from spin correlation functions cal- 
culated by the Lanczos diagonalization method in the 
S* = 1 FS clusters. The first aim is that the LRO of the 
three-sublattice structure is quantitatively confirmed in 
the isotropic case of J2 — Ji by the direct simulation 
of Lanzcos diagonalization in contrast with the situation 
mentioned above, for which no such confirmation in the 
S* = 1/2 case has been successfully made yet. The further 
progress in this study clarifies how the LRO behaves in 
the anisotropic case without having a biased guess based 
on the argument of the classical system. 
Here, we study the Hamiltonian given by 

u^jij2s,-Sj + j2j2 • (1) 

(iJ) [id] 

where Si denotes the S = 1 spin operator. The sum runs 
over all nearest-neighbor pairs having the antiferromag- 
netic interaction of Ji or J2. Energies are measured in 
units of Ji; hereafter, we set Ji = 1. The number of spin 
sites is denoted by Ng. We take Ns/3 to be an integer. 
In all cases of A'g =12, 21, and 27, we have a rhombic 
cluster having an interior angle of 7r/3, at which the two- 
dimensionality may be well captured within FS clusters, 
although only nonrhombic clusters can be formed in the 
cases of Ng —15, 18, and 24 (see Fig. 2). Note that we 
exclude the cases in which a system of three-site chains is 
formed only by Ji bonds at the point of J2 = 0. We im- 
pose the periodic boundary condition for FS clusters. We 
calculate the lowest energy of 7i in the subspace divided 
by J2j = M. The energy is denoted by E{Ng, M). We 
also calculate the correlation function {SfS^) to observe 
the LRO. 

The weak point of the Lanczos diagonalization method 
is that only small clusters can be treated owing to the 




Fig. 2. Shapes of FS clusters of triangular lattice for (a) A^s = 12, 
(b) Ns = 15, (c) Ns = 18, (d) iVa = 21, (e) = 24, and (f) 
TVs = 27. 

exponential increase in the dimension of the Hamiltonian 
matrix with respect to the system size. To overcome this 
problem, we have carried out parallel calculations^^ using 
the MPI-parallelized code^'^ to treat system sizes as large 
as possible. The case of =27 is the largest size in this 
study. In this case, the largest number of dimensions 
is 712 070 156 203 for the subspace of M = 0. Note that 
this number of dimensions is the largest among those 
reported in numerical diagonalization studies of quantum 
lattice systems, to the best of our knowledge. 

Now, we present our results of the S = 1 triangular- 
lattice antiferromagnet for the FS rhombic clusters in- 
cluding the case of — 27; these results are listed in 
Table I. In Fig. 3, the ground-state energy per site Cg — 
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Fig. 3. Ground-state energy per site as a function oi N^' . In the 

1 /2 

inset, the dependence of the spin excitation gap is shown. 

Closed (open) circles represent data for rhombic (nonrhombic) 
clusters. 

E{Ns,M = 0)/Ns is plotted as a function of n!^'^. A 
linear least-squares fitting gives — Cg/Ji = 1.838 ±0.007, 
which is in agreement with the results of the spin- wave 
approximation"^ and coupled cluster method. The in- 
set of Fig. 3 depicts the singlet-triplet energy difference 
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Table I. List of numerical data of the ground-state energy eg = E{Nb,M = 0)/Nb, the singlet-triplet energy difference Ajv^ = 
E{Nb,M = 1) — E(Nb,M = 0), and correlation functions for the rhombic clusters. The distance between sites i and j is denoted 
by I nil- 
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l^il = 2 
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27 


1.870821364074 


0.348963468544 


-0.207869 


0.278656 


-0.125503 


-0.140997 


0.253805 


21 


1.882070035640 


0.441247764571 


-0.209119 


0.286580 


-0.138420 


-0.150456 




12 


1.959106953374 


0.751321679528 


-0.217679 


0.303220 


-0.135128 







A^^ = E{Ns,M = 1) - E{Ns,M = 0). The excitation 
gap seems to vanish in the thermodynamic limit, which 
is consistent with the existence of the LRO. These re- 
sults suggest that the FS clusters shown in Fig. 2 are an 
appropriate series for analyzing our data of correlation 
functions. 

Next, we examine the magnetic order of the triangular- 
lattice antiferromagnet from FS data of not only rhombic 
but also nonrhombic clusters. In the isotropic case, let us 
focus our attention on the signs of {S^Sp in Table I. It 
is possible to group all the spin sites by whether {S^SJ) 
is positive or negative so that sites i and j belong to a 
common group for a positive (SfSj). One finds that the 
grouping divides all the spin sites into three equivalent 
sublattices. Thus, we evaluate the quantity defined as 

a i j^Ai 

where {Si ■ Sj) is evaluated by 3(S'fS'|), some of which 
are obtained from the values of spin correlation func- 
tions presented in Table I, by taking into account the 
isotropic interactions in the spin space. Here, the prime 
at the sum means that j runs over the sublattice Ai in- 
cluding i while excluding the case of j — i. Here, a de- 
notes the label of directions concerned which is chosen as 
Ji among the three directions; we take the average with 
respect to the direction in order to take the nonrhombic 
cases into account. The quantity m^^^^ corresponds to 
the squared sublattice magnetization in the thermody- 
namic limit when the system is magnetically ordered. It 
is known that, in the isotropic case, the spin- wave ap- 
proximation'^ estimates the sublattice magnetization to 
be TO,w = S -Ao + O[l/S] with Aq = 0.261. For S = I, 
TTig^ = 0.546, which will be compared with the value 
extrapolated from our FS w^^^^^ in the isotropic case. 
Note that the quantity nf^^^^ is the same as eq. (2. lid) 
in ref. 7. Note also that the analysis of w^^^^g does not 
directly provide us with information on the relationship 
between spins in different sublattices. 

We depict "^eUag ^^S- 4 as a function of l/Ng^^. 
One can see oscillating behavior in the system size de- 
pendence with respect to whether Ns is odd or even; thus, 
the data series of odd Ng and even Ng should be treated 
separately from each other. Note here that each series of 
'^chag should converge to a common value in the thermo- 
dynamic limit. Then, we perform a linear least-squares 
fitting to each series under the constraint condition of a 
common intercept. From the analysis without Ns = 27, 
we obtain the fitting lines in Fig. 4, drawn as the solid 
line for the data of Ng — 12, 18, and 24 and as the dot- 
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Fig. 4. Extrapolation analysis of FS results for "ij'?j^g in the 
isotropic case. Open diamonds denote the results for iVs = 12, 
18, and 24. Open circles represent the results for Nb = 15 and 
21. The closed circle represents that for A^s = 27, which is the 
largest system in this study. The result of the spin-wave approx- 
imation is plotted by the closed square at the ordinate. See text 
regarding linear lines. 

ted line for the data of iVg = 15 and 21. The common 
intercept is obtained to be 

0.546 ±0.051, (3) 

as an extrapolated value. From the analysis with Ng = 
27, we have the common intercept 

0.551 ±0.035, (4) 

which is in agreement with the result (3). Note that 
both results (3) and (4) are in agreement with the spin- 
wave result mentioned above. Therefore, the constraint- 
condition analysis can provide us with a reliable value 
from extrapolation even without Ng = 27. In Fig. 4, 
the single-dotted and double-dotted chain lines are also 
drawn for even Ng and odd Ng, respectively, as simple 
linear least-squares fittings without considering the con- 
straint condition. One finds that the intercept for even 
Ng underestimates and that for odd Ng overestimates the 
squared sublattice magnetization with which the result 
of the constraint-condition analysis is in agreement. 

Next, we examine the anisotropic case oi J2/J1 < 1; 
the FS result of w^lag ^'^^ each Ng as a function of J2 / Ji is 
shown in Fig. 5. One can observe that the magnetic order 
decreases rapidly as J2 / Ji is decreased. We perform the 
constraint-condition analysis in the cases of J2/J1 = 0.5 
and 0.8; the analysis for extrapolation is shown in Fig. 6. 
The common intercept was obtained to be 0.04 ± 0.17 
for J2/J1 — 0.5, which indicates that the LRO of the 
3SS disappears. On the other hand, we obtain a positive 
common intercept for J2/J1 — 0.8; the result of J2/J1 = 
0.8 is shown in Fig. 5 together with the result in the case 
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Fig. 5. FS results for m^lj^g- Closed squares, diamonds, and cir- 
cles denote the results for A^s = 12, 18, and 24, respectively. 
Open squares and circles represent the results for Nb = 15 and 
21, respectively. Double circles with the error bar denote extrap- 
olated results from the constraint-condition analysis of data up 
to Na = 24 for J2/J1 = 0.8, 0.9, and 1. 
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Fig. 6. Extrapolation analysis of FS results for m^'Jj^ in the 
anisotropic cases of J2/J1 = 0.5 and 0.8. The procedure of the 
extrapolation is the same as that in Fig. 4. 

of J2/J1 = 0.9 and the result (3) of J2/J1 = 1. These 
resuhs suggest that the LRO of the 3SS survives in a 
range having a nonzero width. 

Between J2/J1 = 0.5 and 0.8, unfortunately, the 
constraint-condition analysis gives a quite large error or 
a negative common intercept owing to a severe FS effect; 
it is not so easy to estimate the critial value of J2 / Ji be- 
tween the region where the LRO of the 3SS survives and 
the region where it disappears. In the analysis of only 
even Ng shown in Fig. 7(a), the intercept at the ordinate 
changes its sign between J2/J1 = 0.75 and 0.8. From the 
observation of the isotropic case in Fig. 4, J2/J1 at the 
vanishing intercept is considered to be overestimated as 
the critial value of J2/J1 from the result in Fig. 7(a). On 
the other hand, the analysis of odd Ng shown in Fig. 7(b) 
indicates that the intercept at the ordinate changes its 
sign between J2/J1 = 0.6 and 0.65; J2/J1 at the van- 
ishing intercept is considered to be underestimated as 
the critial value of J2/J1 from the result in Fig. 7(b). 
Therefore, it is reasonable to consider that J2/J1 at the 
boundary is 0.7 ± 0.1. 

In summary, we studied the ground state of an S* = 1 
Heisenberg antiferromagnet on a triangular lattice with 
spatial anisotropy by Lanczos diagonalization calcula- 
tions with large-scale parallelization; the stability of the 



Fig. 7. Extrapolation analysis for m^iag without the constraint 
condition in the anisotropic cases between J2/J1 = 0.8 and 0.5 
with 0.05 intervals for (a) even-Afg results and (b) odd-A^s results. 

long-range magnetic order of the 3SS is examined. We 
found that in the isotropic case, our Lanczos diagonaliza- 
tion data of spin correlation functions are successfully ex- 
trapolated to a long-range-ordered value which is quan- 
titatively consistent with the spin-wave approximation. 
In the present study, we concluded that, in the range of 
0.7 ^ J2/J1 < 1, the sublattice magnetization gradu- 
ally shrinks as J2 is decreased, while the ground state 
maintains the 3SS. The width of this region is narrower 
than those of the magnetically ordered state obtained in 
studies by approximation.^*^' In future studies concern- 
ing this model, the magnetic structure factor as a func- 
tion of wave number should be examined in the range 
near J2/J1 = 1, which will contribute much to our un- 
derstanding of the relationship between the spiral states 
and the states with the LRO of the 3SS obtained in the 
present study. In the region near J2 = at which the Hal- 
dane gap exists, on the other hand, the behavior of the 
spin excitation gap should also be tackled. Since clusters 
up to Ns — 24 take various tilting angles, it is difficult 
to observe systematic behavior of the Ng dependence; 
calculations of the A^s = 27 system with anisotropy are 
required. NiGaS^^ and Ba3NiSb209^° are considered to 
be good candidate materials for the S — 1 Heisenberg 
antiferromagnet on the isotropic triangular lattice, al- 
though experimental observations and theoretical pre- 
dictions are not necessarily in agreement with each other 
in every aspect. Candidate S = 1 materials in anisotropic 
cases would also contribute much to our understanding 
of magnetic phenomena due to frustration. 
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